Fluctuation and dissipation dynamics in fusion reactions 
from stochastic mean-field approach 



o 
o 

(N 
m 



> 
o 

(N 
O 
0^ 
O 



X 



Sakir Ayik,^ Kouhei Washiyama,^ and Denis Lacroix^ 

Physics Department, Tennessee Technological University, Cookeville, TN 38505, USA 
^GANIL, Bd Henri Becquerel, BP 55027, U076 Caen Cedex 5, France 
(Dated: February 3, 2009) 

By projecting the stociiastic mean-field dynamics on a suitable collective path during the entrance 
channel of heavy-ion collisions, expressions for transport coefficients associated with relative distance 
are extracted. These transport coefficients, which have similar forms to those familiar from nucleon 
exchange model, are evaluated by carrying out TDHF simulations. The calculations provide an 
accurate description of the magnitude and form factor of transport coefficients associated with 
one-body dissipation and fluctuation mechanism. 
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I. INTRODUCTION 

The self-consistent mean-field theory, also known as 
time-dependent Hartree-Fock (TDHF), by employing 
Skyrme-type effective interactions, has been extensively 
applied to describe nuclear collision dynamics at low 
bombarding energies below 10 MeV per nucleon P, 0, 
H, In the mean-field theory, short range two-body 
correlations are neglected and nucleons move in the self- 
consistent potential produced by all other nucleons. This 
is a good approximation at low energies since Pauli block- 
ing is very effective for scattering into unoccupied states. 
Consequently, in the mean-field theory, collective energy 
is converted into intrinsic degrees of freedom via interac- 
tion of nucleons with the self-consistent mean-field, so- 
called one-body dissipation [5i, ^Od . One-body dissipation 
mechanism plays dominant role in low energy nuclear dy- 
namics including deep-inelastic heavy-ion collisions and 
heavy-ion fusion reactions. One important limitation of 
the mean-field theory is related with dynamical fluctu- 
ations of collective motion. In the mean-field descrip- 
tion, while single-particle motion is treated in quantal 
framework, collective motion is treated almost in clas- 
sical approximation. Therefore, TDHF provides a good 
description for average evolution, however it severely un- 
derestimates fluctuations of collective variables. 

On the other hand, it is well known that no dissipa- 
tion takes place without fluctuations 0, [l] . Much effort 
has been done to improve one-body transport descrip- 
tion beyond the mean-field. Most of these transport de- 
scriptions take into account dissipation and fluctuation 
mechanisms due to two-body collisions, which play an 
important role in nuclear dynamics at intermediate en- 
ergies d, [13, ■ Here, we deal with nuclear dynam- 
ics at low energies at which one-body dissipation and 
associated mean-field fluctuations play a dominant role 
in dynamical evolution. One of the fundamental ques- 
tions is how to improve the mean-field theory by incor- 
porating one-body fluctuation mechanism at a micro- 
scopic level? In a recent work, based on an appealing 
idea of Dasso [H, , this question has been addressed. 



A stochastic mean-field (SMF) approach has been pro- 
posed for describing fluctuation dynamics [l^, [l^ • For 
small amplitude fluctuations, this model gives a result for 
dispersion of a one-body observable that is identical to 
the result obtained through a variational approach [T7| . 
It is also shown that, when the SMF is projected on a 
collective variable, it gives rise to a generalized Langevin 
equation [Tsj . which incorporates one-body dissipation 
and one-body fluctuation mechanisms in accordance with 
quantal dissipation-fluctuation relation. These illustra- 
tions give a strong support that the SMF approach pro- 
vides a consistent microscopic description for dynamics 
of density fluctuations in low energy nuclear reactions. 
In this paper, we present another demonstration of the 
SMF approach. 

In a recent work, by a suitable definition of collective 
variables of relative motion, the nucleus- nucleus potential 
energy and one-body friction coefficient as a function of 
relative distance have been extracted from simulations of 
microscopic TDHF [T^, see also [2^. Such a reduction 
is not constrained by adiabatic or diabatic approxima- 
tion, therefore it should provide an accurate description 
of conservative nucleus-nucleus potential energy and the 
magnitude of the one-body dissipation mechanism [2l| . 
It is of great interest to deduce magnitude of diffusion 
coefficients associated with collective variables. How- 
ever, this information is not contained in the standard 
mean-field approximation. The SMF approach provides 
a proper framework for extracting dissipation and fiuc- 
tuation properties of collective variables. In this work, 
we carry out a similar macroscopic reduction of the SMF 
approach on a collective path. In this manner, we deduce 
not only nucleus- nucleus potential and one-body friction 
coefficient, but also one-body diffusion coefficients asso- 
ciated with collective variables. 

In Sec. II, we give a brief description of the SMF ap- 
proach. In Sec. Ill, we present a suitable definition of 
collective variables in heavy-ion collisions, and the cor- 
relation function of Wigner distribution. In Sec. IV, we 
derive transport coefficients associated with relative mo- 
tion from the SMF approach. In Sec. V, conclusions are 



2 



given. 



II. STOCHASTIC MEAN-FIELD APPROACH 

In the standard TDHF, temporal evolution of the sys- 
tem is described by a single Slater determinant con- 
structed with time-dependent single-particle wave func- 
tions ^jaT{r,t). Evolution of single-particle wave func- 
tions are determined by the TDHF equations with proper 
initial conditions, 

d 

ih—<^j„r{r,t) = h{p)<^j„r{r,t), (1) 

where h{p) denotes the self-consistent mean- field Hamil- 
tonian with the one-body density p. For clarity of 
notation spin-isospin quantum numbers t ~ (proton, 
neutron) and a — (spin-up, spin-down) are explicitly in- 
dicated in these expressions. In many situations, it is 
more appropriate to express the mean-field approxima- 
tion in terms of the single-particle density matrix, 

p(r,r',i) (2) 

jar 

where nj'^ denotes occupation factors of single-particle 
states. In the standard TDHF, occupation factors are 
one and zero for the occupied and unoccupied states, 
respectively. If the initial state is at a finite tempera- 
ture, the average occupation factors are determined by 
the Fermi-Dirac distribution. 

TDHF provides a deterministic evolution of the single- 
particle density matrix, starting from a well-defined ini- 
tial state and leading to a well-defined final state. In 
order to incorporate fluctuation mechanism into dynam- 
ics, we give up standard description in terms of a single 
Slater determinant, and consider a superposition of de- 
terminantal wave functions. As a result of correlations, 
initial density cannot have a deterministic shape, but it 
must exhibit quantal zero-point fluctuations, and if the 
initial state is at a finite temperature, it also involves 
thermal fluctuations. In the SMF approach the initial 
density fluctuations are incorporated into the description 
in a stochastic manner [l^. The initial density fluctua- 
tions are simulated by representing the initial state in 
terms of a suitable ensemble. In this manner, an ensem- 
ble of density matrices is generated, 

/(r,r',t) = ^$*^,(r,t;A)4(ar)<I>,,,(r',i;A). (3) 

ijar 

Here '^jcrri'i'.t] A) is a complete set of single particle ba- 
sis, A denotes a member in the ensemble, and matrix 
elements p^j (crr) are time-independent random Gaussian 
numbers. Gaussian distribution of each matrix element 
is specified by a mean value p^- (ctt) = Sijn'j'^ , and a vari- 



ance, 

p'M)p^'A<^'r') 

= ls,,,6u'Srr'S„„, [nr(l - nf ) + (1 - nD] , 

(4) 

where n°"^ denotes the average occupation factor for 
a given values of spin-isospin quantum number ct and 
T. Stt' a-nd <5ctct' indicate that density matrix elements 
are assumed to be uncorrelated in spin-isospin space. 
A member of the ensemble is generated by evolving 
the single-particle wave functions according to the self- 
consistent mean-field of that member, 

d 

ih—^j,r{r,t;\) = h{p^)<^,,r{r,t-\), (5) 

where h{p^) is the self-consistent mean- field Hamiltonian 
in that event. 

HI. STOCHASTIC WIGNER DISTRIBUTION 

In order to carry out projection of the SMF on a collec- 
tive space, to determine transport coefficients of collec- 
tive variables and to establish connection with the collec- 
tive transport models, it is very convenient to introduce 
the stochastic Wigner distribution. The Wigncr distri- 
bution for each event A is defined as a partial Fourier 
transform of density matrix as 

X E "^J- + f ' ^' ^) PU'^^^^^-r (r - |, t; a)(.6) 

ijtTT 

In this work, we focus on head-on collisions of two 
heavy-ions and take the collision direction as the a;- 
axis. Following Ref. [l^, we define center-of-mass coor- 
dinate , total momentum and mass number Aj_ of 
projectile-like (-I-) and target-like (— ) fragments by intro- 
ducing the separation plane. The separation plane can be 
conveniently defined as the plane at position where iso- 
contours of projectile-like and target-like densities cross 
each other. We indicate position of the separation plane, 
i.e., position of the window a,t x = xq. Illustration of den- 
sity profiles and separation plane locations are displayed 
at different times of the symmetric reaction ''"Ca-l-^^Ca 
in Fig. [TJ For calculations in this figure and the rest 
of the paper, we use three-dimensional TDHF code de- 
veloped by P. Bonche and co-workers with the SLy4d 
Skyrme effective force [l^ and for technical details please 
see Ref. 

It is convenient to express macroscopic variables in 
each event in terms of the reduced Wigner distribution 
f^{x,px,t) according to 

Ri = ^j ^Oixo±x)xf\x,px,t), (7) 
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FIG. 1: Density profiles p{x,y,0) obtained with TDHF for 
the ''°Ca-l-''°Ca reaction at Sc.m. = 100 MeV at different R. 
The iso-densities are plotted every 0.025 fm~"^. In each case, 
the vertical line indicates the separation plane. 



and 



dxdpi 

2TTh 



dxdpx 

2TTh 



e{xo±x)f^{x,p,,t). 



(8) 



(9) 



We note that these definitions do not involve semi- 
classical approximations and are fully equivalent to those 



given m 



The ratio P^/R 



(R) determines 

inertia of both sides of the window and the relative mo- 
mentum is defined as 



where ii^{R) 



^A pA 



H^{R)R^ 



lm^/{m^ + m^) and R' 



(10) 



are the reduced mass and the relative velocity of projec- 
tile and target sides, respectively. The reduced Wigner 
distribution f^{x,Px,t) is obtained by integrating over 
the phase-space variables y, z,Py,Pz according to 



f\x,p^,t) 



(11) 



In order to extract diffusion coefficients associated with 
collective variables, we need diffcrcnt-time correlation 
function of the reduced Wigner distribution on the win- 
dow. Assuming that the amplitude of density ffuctua- 
tions is small, this correlation function on the window 
is calculated in the semi-classical approximation in Ap- 
pendix A to give 



Sf^{x,Px,t)df^{x,p'^,t')\^^^„ 

TTl 

{2TTh)—6{px-p'^)6{t-t')A+{xo,Px,t), (12) 

\Px\ 



where the quantity {xq ,Px,t) is defined as 
A^ixo,Px,t) = ^{/^^(xo,Px,i) [l-/r(a^o,Px,i)] 

±f^^Xo,Px,t) [l-f^^Xo,Px,t)]}. 

(13) 

In this expression, fP'{x,Px,t) denotes, in spin-isospin 
channel (cr, r) , the average value of reduced Wigner func- 
tion associated with wave functions originating from pro- 
jectile, 

X 51*"^^ (^+ Y'?/'^'^) ^r^iar {x-^,y,z,t 



The average quantity 

ff.^{xo,Px,t) 



ff.^xo,px,t) 
n{xo,t) 



(14) 



(15) 



denotes the reduced Wigner distribution averaged over 
phase-space on the window, i.e., on the plane divid- 
ing projectile-like and target-like nuclei, where ^l{xo,t) 
is the phase-space volume over the window. Quantities 
It^ {xo,Px,t) and fx^ {xo,Px,t) are average values of re- 
duced Wigner function associated with wave functions 
originating from target in spin-isospin channel (cr, r) , 
which are defined in a similar manner. 

We approximate the phase-space volume over the win- 
dow as 



n{xQ,t) = Trrl^^i^{xo,t) 



(16) 



where pp stands for the Fermi momentum. In this ex- 
pression rneck{xQ,i) dcuotes the equivalent sharp radius 
of the neck, which is defined as 



neck 



{xo,t) 



1 



no(xo,t) 



dydzn{xo,y, z,t), (17) 



where n(xo,y, z,t) is the local nucleon density while 
no{xo,t) denotes the density at the center of the neck, 
i.e., no(xo,t) = n(xo,0, 0,t). The evolution of de- 
duced from Eq. p?)) is shown by solid lines in Fig. |2] for 
the '^OCa+^OCa reaction as a function of relative distance. 
While the neck radius has rather reasonable values at 
small R, Eq. P7)l leads to unrealistic large values for well 
separated nuclei. To overcome this difficulty, we use an 
alternative approach by considering that fp'^rp{xQ,Px,t) 
should be close to 1.0 around the average value of p^;. 
By imposing this condition, we directly determine an 
approximate phase-space volume from Eq. IjlSp . Then, 
we deduce r„eck at each relative distance R by inverting 
Eq. ((TBI) . These are indicated by filled circles in Fig. |21 
We see that the second prescription not only provides 
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FIG. 2; Neck radius determined for the '*''Ca+*'^Ca reac- 
tion at Ec.m. = 100 MeV by using Eq. (fT7)l (solid line), and 
by imposing the condition that the reduced Wigner function 
/p/y(a;,pa;, t) is close to 1.0 around the average momentum 
(filled circles). 



a reasonable behavior for Tncck at large distances, but 
also matches Tnock deduced by using Eq. pT)) at small 
distances. In the calculations we use the effective neck 
radius determined by the second approach. Examples of 
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FIG. 3: Reduced Wigner function p" {xo,px) averaged over 
phase space on the window [x = a;o) for projectile-like (solid 
line) and target-like (dotted line) nuclei for the '"^Ca-|-'*''Ca 
reaction at -Ec.m. = 100 MeV at different R. The Fermi mo- 
mentum pf is taken as 270 MeV/c. 



reduced Wigner function are shown in Fig. [3] for different 
relative distances. Not surprisingly, the reduced Wigner 
function is sometimes above 1 or below 0. This is in- 
deed expected since the full quantum Wigner transform 
is considered without making use of any semiclassical ap- 
proximation. 



IV. MOMENTUM DIFFUSION COEFFICIENT 

In a recent work ^19J, considering simple one- 
dimensional macroscopic reduction of TDHF, average 
transport properties of relative motion in heavy-ion col- 
lisions have been investigated. Temporal evolution of av- 
erage value of relative distance R{t) = R^{t) and aver- 
age value of relative momentum P{t) = P^{t) are calcu- 
lated for average mean-field trajectory determined by the 
TDHF wave functions. Relative motion of colliding ions 
were analyzed in the basis of a simple classical equation 
of motion, 



(18) 



Knowing time evolution of R{t) and P{t), average col- 
lective properties, namely, nucleus-nucleus potential en- 
ergy U{R) and form factor of one-body friction coefficient 
7(-R) are determined by inverting Eq. psp . In this work, 
we consider the same geometry of head-on collision of 
heavy-ions and extend the macroscopic reduction treat- 
ment by considering the SMF approach. We analyze the 
relative motion by employing a Langevin equation. The 
Langevin equation for the relative motion has the form. 



dt 



dR^ 



U{R')-jiR^)R'+ep{t), (19) 



where (t) is a Gaussian random force acting on the rel- 
ative motion. Ignoring non-Markovian effects, the ran- 
dom force reduces to white noise specified by a correlation 
function, 



2S{t-t')Dpp{R). 



(20) 



Here Dpp{R) denotes the momentum diffusion coeffi- 
cient, which may depend on the mean value of the relative 
distance R. 

In order to extract momentum diffusion coefficient, we 
calculate the rate of change of the relative momentum 
employing the SMF equations. The rate of change of 
relative momentum involves kinetic parts due to nucleon 
exchange between projectile and target, and also involves 
terms arising from potential energy. In the previous in- 
vestigation [2l[ , it is observed that during evolution from 
the entrance channel until passing over the Coulomb bar- 
rier, one-body dissipation mechanism is strongly corre- 
lated with nucleon exchange between projectile-like and 
target-like nuclei. This behavior is similar to phenomeno- 
logical nucleon exchan ge m odel and the window formula 
for energy dissipation [23l. [2^ . Therefore, in the equa- 
tion for the rate of change of relative momentum, we 
consider only kinetic terms corresponding to momentum 
flow across the window, which can be conveniently ex- 
pressed in terms of reduced Wigner distribution as 



dt 



-P^ 



dpx pI 
2'Kh m 



f^{x,px,t)\x=xo + potential terms. 



(21) 
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Small fluctuations of relative momentum are connected 
to small amplitude fluctuations in Wigner distribution. 
Ignoring contribution arising from potential terms, we 
have for small fluctuations of relative momentum 



dt 



-SP' 



dpx pI 
2'Kh m 



epit)- (22) 



The right hand side in this expression acts as a random 
force for generating fluctuations in the relative momen- 
tum. Since Sf^{x,px,t) is a Gaussian random quantity, 
the random force ^p(t) is also Gaussian random, which 
is specified by a correlation function, 



dpx dp'^ pI p'l 



2'Kh 2nh m m 



x5f^{x,px,t)5f\x,p'^,t')\,=x,m) 

Using the expression for the correlation function of the 
reduced Wigner distribution in Eq. \12\ . according to 
Eq. (|20p. the momentum diffusion coefficient is given by 



Dpp{t) 



dpx \Px\ Px A+/ ,\ ^r,A\ 

TT^ -TT^ i^O, Px,t). (24) 

ZTrh m Z 



From the SMF approach, we cannot directly derive an 
expression for the friction coefficient j{R). The reason is 
that we cannot associate the net momentum flow across 
the window, which is given by the first term on the right 
side of Eq. (PT|) , with dissipative force acting on the rel- 
ative motion. However, from the expression (|24|) for dif- 
fusion coefficient and from the random walk mechanism 
of nucleon exchange [2^, !2i| , we can infer an expression 
for the friction coefficient. In the expression for diffusion 
coefficient, first and second terms correspond to nucleon 
flux from projectile to target and from target to pro- 
jectile, respectively. Each nucleon transfer changes the 
relative momentum by an amount p,j; and increases the 
dispersion of the relative momentum by an amount p^. 
Nucleon transfer in both direction increases dispersion 
of relative momentum, therefore diffusion coefficient is 
determined by total nucleon flux, i.e., sum of flux from 
projectile to target and from target to projectile. On 
the other hand, dissipation is determined by the net mo- 
mentum flow through the window. Hence, the resultant 
dissipative force can be expressed as 



F{t)^ I ^^P^PxA-ixo,Px,t). 



(25) 



Then, it is possible to deduce from TDHF simulations the 
momentum diffusion coefficient Dpp{t) = Dpp{R) and 
the friction force F{t) as a function of relative distance. 
We note that these transport coefficients correspond to 
the phenomenological window formula arising from the 
nucleon exchange mechanism j23], and they are deter- 
mined in terms of the average evolution specified by the 
TDHF. 

Rather than calculating the dissipative force, it is more 
instructive to calculate the friction coefficient j{R)- For 



this purpose, we assume that dissipative force is propor- 
tional to relative velocity, i.e., F{t) = —^{R)R, and con- 
sider the reduced friction coefficient (3{R) = j{R)/iJ,{R), 
where fJ-{R) denotes inertia associated with relative mo- 
tion. Solid lines in Fig. |4] show the reduced friction 
coefficient as a function of R for head-on collision of 
^"Ca-l-^^Ga at two different center-of-mass energies. For 
each energy, enlarged plot around the Coulomb barrier 
region is shown in the insert. In a recent work [2]| . 
we extracted the reduced friction coefficient associated 
with relative motion employing a different reduction 
procedure, so-called Dissipative-Dynamics TDHF (DD- 
TDHF), which, in principle, incorporates dissipation due 
to both window and wall mechanisms. Dashed lines in 
Fig.[5]show the results of this reduction procedure. Good 
agreement is found between two different calculations 
above and close to the Coulomb barrier (~ 9.8 fm). Be- 
low the Coulomb barrier, the DD-TDHF method is not 
reliable. However, the method based on the SMF pro- 
vides a proper description of the one-body friction co- 
efficient due to nucleon exchange mechanism for a wide 
range of relative distance. 
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FIG. 4: Reduced friction coefficient /3(_R) = j{R)/iJ, as a 
function of R for the "^^Ca-t-^^Ca reaction at Bc.m. = 100 MeV 
(upper panel) and at Ec.m. ~ 150 MeV (lower panel). For 
each energy, a zoom on the Coulomb barrier region is also 
shown in the insert. 



Solid lines in Fig. [5] show the momentum diffusion co- 
efficient Dpp, Eq. ([M)) . as a function of R for head-on 



collision of *°Ca+^°Ca at two different center-of-mass en- 
ergies. Similarly to the reduced friction coefficient, mag- 
nitude of the momentum diffusion coefficient increases for 
decreasing relative distance. The increase of magnitude 
of transport coefficients, i.e., friction and diffusion coef- 
ficients, for decreasing R is essentially due to larger win- 
dow area and larger number of nucleon exchange between 
projectile-like and target-like nuclei. It is important to 
realize that, even though the ordinary TDHF does not 
contain information about density fluctuations, we can 
employ the average information provided by the TDHF 
to calculate diffusion coefflcients associated with macro- 
scopic variables. In practical applications, the momen- 
tum diffusion coefflcient is usually taken as the thermal 
equilibrium value determined by the Einstein relation in 
terms of friction coefficient and effective temperature as 



: ^{R)T{t) 



(26) 



In this expression, T{t) denotes the effective tempera- 
ture assuming local equilibrium. It can be determined 
in terms of excitation energy denoted by E* by the re- 
lation T{t) — ^J E*{t) / a, where a denotes level density 
parameter, taken here as a = A/12. We can estimate the 
excitation energy in terms of dissipated energy according 
to 



Jo 



(27) 



Dashed lines in Fig. [5] show the diffusion coefficient 
D'Pp(R) determined according to the Einstein relation. 
As seen from the figure, the Einstein relation severely 
underestimates magnitude of dynamical diffusion coeffi- 
cient. The fact that the Einstein relation severely under- 
estimates the dynamical diffusion coefficient associated 
with the relative motion was already realized in the phe- 
nomenological description of nucleon exchange model in 
Ref. M- 



V. CONCLUSIONS 

Recently proposed stochastic mean-field theory incor- 
porates both one-body dissipation and fluctuation mech- 
anisms in a manner consistent with quantal fluctuation- 
dissipation theorem of non-equilibrium statistical me- 
chanics [l^. This was illustrated for slow collective mo- 
tion by projecting equation of motion of the SMF onto 
a collective space in adiabatic limit. The projection 
gives rise to a generalized Langevin equation for collec- 
tive variables, in which mean-fleld dissipation and fluc- 
tuation mechanisms are connected through the quantal 
fluctuation-dissipation relation. Therefore, this approach 
provides a powerful tool for microscopic description of 
low energy nuclear processes in which two-body dissipa- 
tion and fluctuation mechanisms do not play important 
role. The low energy processes include induced fission. 
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FIG. 5: Diffusion coefficient Dpp obtained by SMF (solid 
lines) and by the Einstein relation D'pp — 'y{R)T{t) (dashed 
lines) as a function of R for the "Ca-l-^'^Ca reaction at 
-Ec.m. = 100 MeV (upper panel) and at Ec.m. = 150 MeV 
(lower panel). 



heavy- ion fusion near barrier energies, and spinodal de- 
composition during the expansion phase of hot piece of 
nuclear matter produced in heavy- ion collisions [16. . ,25] . 

In this work, we carry out a different projection of 
the SMF approach on the relative motion in fusion re- 
action by following the DD-TDHF method introduced 
in [l9| and deduce one-body friction and one-body dif- 
fusion coefficients associated with relative motion. It 
is remarkable that expressions of transport coefficients 
for the relative motion (as well as transport coefficients 
for other macroscopic variables which are not mentioned 
in this work) have the same form as given b y th e phe- 
nomenological nucleon exchange model [1^ |24| . The 
phenomenological nucleon exchange model involves an 
important assumption, namely, when a nucleon passes 
through the window, it instantaneously equilibrates with 
the new environment on the other side of the window. 
On the other hand, transport coefficients deduced from 
the SMF approach do not involves this assumption, and 
also they are not restricted by adiabatic or diabatic ap- 
proximation. Therefore, these transport coefficients pro- 
vide a microscopic basis for determining magnitude of 
the actual one-body dissipation and the corresponding 
mean-field fiuctuation mechanism. We also stress the 
fact that, assuming amplitude of density fluctuations are 
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small, transport coefficients are calculated in terms of the 
average evolution determined by TDHF simulations as a 
function of relative distance. In the continuation of the 
investigations, we plan to generalize the projection pro- 
cedure of the SMF approach for off-central collisions and 
also deduce transport coefficients for nuclcon diffusion in 
grazing heavy-ion collisions. 
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APPENDIX: CORRELATION FUNCTION OF 
WIGNER DISTRIBUTION 

Small amplitude fluctuations of Wigner distribution 
can be expressed as 

ij(TT 

(A.l) 

I 



where the single-particle wave functions are complete set 
of solutions of the ordinary TDHF. The initial values 
of stochastic expansion coefficients Spfj{aT) are Gaus- 
sian random numbers as specified by Eq. In princi- 
ple, these coefficients evolve in time according to time- 
dependent RPA equations. Here, we ignore this evo- 
lution and take them as Gaussian random numbers as 
specified by the initial conditions. Fluctuating part of 
the density matrix can be separated into four groups, 
Sppp, Sp^rp^ SppT ^-iid Sp^p, which are associated with 
wave functions originating from projectile and target nu- 
clei and the mixed terms. As a result, small ampli- 
tude fluctuations of the Wigner distribution separate into 
four parts, Sf^p{r,p,t), 5f^j.{r,p,t), Sf^j,{r,p,t) and 
6f^p{r,p,t). We calculate the equal time correlation 
function of the Wigner distribution in semi-classical ap- 
proximation. First, we consider the correlation function 
associated with wave functions originating from projec- 
tile. Using the expression ^ for the variance of the 
matrix elements, we deduce 



Sf^p{r,p,t)Sf^p{r',p',t) 



d^sd^s' exp ( ~^p ■ s I exp ( — — p' • s 



(A.2) 



In the term that is proportional to we use the closure relations to find, 



i£P 



S + S 



(A.3) 



In this expression, summation runs over a complete set of single-particle states, i.e., occupied and unoccupied states 
originating from the projectile. The closure relation satisfied by the complete set of states at the initial state will 
remain valid at later times. In the second step, we introduce the Wigner distribution. 



Y.^U (r + -,t)nj-$,..(r'--,t) = 



■ exp 



s + s 



Q 



/err 

IP 



r + r s — s 



Q,t], 



(A.4) 



where fp'^{r,p, t) denotes the ensemble averaged Wigner distribution associated with wave functions originating from 
projectile in spin-isospin channel (cr, r). After making transformations, s = -I- 77 -I- Y/2 and s' — — ij + Y/2, the 



8 



term that is proportional to occupation factor nj"^ in the right-hand side of Eq. (|A.2p becomes 



J2 d^Yd^T^exp 



P + P 



exp 



^-T]- {p-p') 



Y ri 



fp [r-^ + ^.Q,t]5[r' -r + - 



(A.5) 



Assuming that the Wigner distribution is a smooth function of r, — T + i ' ~ /p^(''' ^) ^'^'^ ^('''~'''+t) ~ 
— r'), we can carry out the integrations over r] and Y to obtain (27r/i)'^(5(p — p') and {2'KKf' 5{p — Q) , respectively. 
As a result, the term ()A.5|1 becomes 



(A.5) - (2^;i)35(p - pO-^lr - r') ^ /i^-(r, p, t). 



(A.6) 



For the term proportional to nf'^n'j'^ in Eq. (|A.2p . again we introduce the Wigner distribution for the factor involving 
the index j, 



jeP 



and for the one involving the index i, 



■ exp 



s + s 



Qi 



far 
JP 



r + r s — s 



(A.7) 



r — r + 



Q2 



IP 



r + r s — s 



-,Q2,t 



(A., 



Making the same transformations, s = +77 + Y/2 and s' ~ t] + Y /2, the term that is proportional to nf'^n'j'^ in the 
right-hand side of Eq. (jA.2[) becomes 



E 



X exp 



d-^Yd'^rjexp 
Y 



r r 



Q2 



Y 



P + F 



exp 



-11 {P- P') 



d^Qi d^Q2 

{2TThf {2lThY 



exp 



Qi 



r + r T7 



2 -2'^''* 



(A.9) 



We introduce another change of variables Qi = Q + q/2 and Q2 = Q q/2, and again assume that the Wigner 
distribution has a smooth function of r and ignore r] dependence. Then, integrations over 77, q and Y give {2Trh)^S{p — 
p'), {2TTh)^6{r — r') with Qi « Q2 = Q and {2Trh)^d{p — Q), respectively. As a result, the term (|A.9p becomes 



(A.9) = (2^fi)=^5(p-p')'5(r-r')E/p'(^'P'*)/p'(^'P'*)- 



(A.IO) 



Combining together, equal time correlation function (IA.2|) of the Wigner distribution associated with wave functions 
originating from projectile becomes. 



Sf^pir,p, t)Sf^p{r',p', t) = {2nhf5{p - p')5{r - r') ^ fp^r^p, t)[l - f^^r,p, t)]. 



(A.ll) 



In a similar manner, we can calculate the correlation function of the Wigner distribution associated with wave functions 
originating from target and from mixed configuration. 



6f^^{r, p, t)Sf^^ir',p', t) = {2nhf5{p - p')5{r - r') ^ /?^(r,p, t) [1 - /?-(r, p, t)\ (A.12) 



and 



6f^^{r,p,t)5f^T{r',p',t) = {2^hf6{p - p')5{r - r')K+{r,p,t), 



(A.13) 
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where 

A+(r,p,t) = Y,{fp{r,p,t) [1 - f^^{r,p,t)]+f^^r,p,t) [1 - f^^r,p,t)]} . (A.14) 

err 

Total correlation function of the Wigner distribution is the sum of (jA.ll[) . (|A.12[) and (|A.13[) . In the mean- 
field description, the sub-spaces of wave functions originating from projectile and target nuclei behave like pure 
states. Therefore, contributions of correlations coming from direct terms involving /p^(t,P, [1 — fp'{i',p,t)] and 
f^'^ {r,p,t) [1 — fj''^{r,p,t)] are expected to be small. Hence, we can approximately express the total correlation 
function of Wigner distribution as, 

Sf^{r,p,t)6f^{r',p',t) « {2T:hfS{p - p')d{r - r')A+{r,p,t). (A.15) 

We also want to calculate different time correlation function of the Wigner distribution. Assuming that the correla- 
tion function has short correlation time, i.e., much shorter than mean- free path, different time correlation function can 
be deduced by observing that in short time intervals of order of correlation time |i — i'| < Tcorr, Wigner distribution 
may be approximated as a free propagation, 6f{r,p, i + r) « Sf{i — Tp/m,p, t). As a result, different time correlation 
function can be expressed as, 

6f^{r,p,t)Sf^{r',p',t') = {2nhf5{p~p')S[r-r' - {t - t')p/m]A+{r,p,t). (A.16) 

In order to deduce the correlation function on the window, x — x' = xq, we notice that 

6[r-r' -{t- t')p/m] ^ —S{t - t')S{y - y')S{z - z'). (A.17) 

\Px\ 

In determining transport coefficients, we need to carry out integration over window variables, y, z^Py,Pz, of product 
of Wigner distributions. Since construction of three-dimensional Wigner functions in terms of TDHF wave functions 
requires a large numerical effort, we introduce the following approximation for the phase-space integration over the 
window, 

/ [ dydz'^^^rp^{r,p,t)f^^{r,p,t) « (A.18) 



Here Vl{x^t) denotes the phase-space volume on the win- can be expressed in terms of the reduced Wigner distri- 
dow. As a result, the correlation function on the window butions along x-axis given by Eq. 
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